Interpretation of the de Haas - van Alphen experiments in MgB 2 



I. I. Mazhf] 

Code 6391, Naval Research Laboratory, Washington, DC 20375-5000, USA 

Jens Kortus^ 

Max- Planck- Institut fur Festkorperforschung, Heisenbergstr. 1, D-70569 Stuttgart, Germany 

(Dated: February 1, 2008) 

Recent reports on quantum oscillations in MgB2 provide valuable information on three important 
aspects of this material: (1) electronic structure near the Fermi level, (2) disparity of the electron- 
phonon interaction between the two systems of bands (3) renormalization of spin susceptibility. 
However, extraction of most of this information requires highly accurate band structure calculations 
of the relevant quantities. In this paper we provide such calculations and use them to analyze the 
experimental data. 



MgB2, a novel superconductor with T c w 40 K, has 
attracted enormous attention in the last-,year. The most 
popular model, suggested by Liur-,et aln and Shulga el 
al.p and elaborated by Choi et al.jfl is the two-gap model, 
which, based on the very large interband disparity of the 
electron-phonon interaction (first noted in Ref.^), pre- 
dicts two different ganSj-for the two different band sys- 
tems. The calculationsotj yield an effective (including an 
enhancement due to gap variation) electron-phonon cou- 
pling constant of the order of 1. On the other hand, the 
two-gap theory has a serious conceptual problem: Two 
distinctive gaps may exist only if the interband impurity 
scattering is very weak. That seems to be in contrast 
to the experimental observation that even poor quality, 
high-resistivity, samples have vera good superconducting 
properties. It has been arguedfm that the specificity of 
the electronic and crystal structure of MgB2 results in 
a peculiar relation among the three relevant relaxation 
rates, namely that the impurity scattering inside the so- 
called 7T band is much stronger than inside the a band, 
and the latter, in turn, is much stronger than the in- 
terband scattering. However, there has been no direct 
experimental confirmation of this claim. 

On the other hand, some authorso argue that the calcu- 
lated band structure is strongly renormalized by electron- 
electron interactions not accounted for in the local den- 
sity calculations, so that the plasma frequency is a factor 
of five smaller than the calculated one. This would im- 
ply an electron-phonon coupling constant less than 0.2. 
There are claiais that infrared spectroscopy supports this 
point of view Jala although other researchers in the fieldEl 
dispute the interpretation accepted in Refs.||,^[ In any 
case, the fact that all optical experiments till now have 
been performed on polycrystalline samples, undermines 
their value as a decisive test for the electronic structure 
calculations. 

The first single crystal angular-resolved photoemis- 
sion (ARPEJI)i9 measurements agree very well with the 
calculations xA However, some calculated bands have not 
been observed, and, furthermore, ARPES probes only a 
very thin surface layer and is therefore often not repre- 
sentative of the bulk electronic structure. 

Historically, the most reliable probe of the bulk elec- 



tronic structure has been the de Haas-van Alphen effect 
(dHvA)_Recent observation of this effect in MgB-2 single 
crystals^ provides key information to assess the validity 
of the standard band structure calculation. Given the 
fact that most theoretical papers rely on this band struc- 
ture, the importance of a proper analysis of these data 
can hardly be overestimated. It must be emphasized that 
such an analysis requires highly accurate band structure 
calculations, i.e., the use of a much finer k-point mesh in 
the Brillouin zone and a much more accurate integration 
than is customary in other applications of the band the- 
ory. In this paper we present such calculations and show 
that both Fermiology and effective masses (and hence 
the Fermi velocities and plasma frequencies) produced 
by conventional band structure calculations are in excel- 
lent agreement with the experiment, thus giving a strong 
foundation for the widespread use of this band structure. 
Furthermore, we show that the calculational predictions 
of a strong disparity of the electron-phonon interaction 
in the two band systems in MgB2 are supported by the 
de Haas- van Alphen experiment, and that the scattering 
rates inside the a band and between a and 7r bands are 
probably much smaller than inside the 7r bands. 

The Fermi surface of MgB2 consists of four sheets.E^I 
Two sheets come primarily from the boron p x and p y 
states, and form slightly (nearly sinusoidalUd warped 
cylinders, a (bonding) and a* (antibonding) ,L3 and two 
tubular networks, the bonding one, n, in the T (k z = 0) 
plane, and the antibonding one, 7r*, in the A (k z = tt/c) 
plane. There are 6 extremal cross-sections for the field 
parallel to k z (along the T— A line). These are: (1) a in 
the r plane; (2) a* in the T plane; (3) ir in the T plane 
("holes" between the tubes); (4) a in the A plane; (5) a* 
in the A plane; and (6) ir* in the T plane. For a field 
parallel to k y (perpendicular to the T— AM plane) there 
are two extremal cross-sections (tubes' necks), for the it 
surface (7) and for the ir* surface (8). 

We performed highly accurate and well converged full 
potential linear augmented plane wa^je, (LAPW) calcu- 
lations, i-tising the WIEN-97 package,t3 including local 
orbitalsli3 to relax the linearization errors. We used the 
Generalized Gradient Approximation of Perdew-WangEH 
for the exchange-correlation potential. By comparing the 
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results with LMTO (linear muffin tin orbitals) calcula- 
tions, we found that for a proper description of the a 
orbits it is essential to use a full potential method. It is 
furthermore essential to use a very fine mesh in k-space; 
we employed a 38x38x27 mesh, corresponding to 1995 
inequivalcnt k-points. To achieve sufficient accuracy for 
the small areas of the orbits 1,2,4,5,7 and 8, we used, 
an integration engine built in the SURFER program,E3 
which internally interpolates the integrand with splines. 

The bare (band) masses in the third column of table 
I were then calculated by varying the Fermi energy and 
using the standard formula, 
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Here and below we use the notation A for the areas of the 
orbits in standard units and F for those in Tesla units. In 
order to obtain the energy derivatives we fitted the cal- 
culated A(E) by quadratic polynomials in the ranges of 
about 0.03 Ry around the Fermi energy. The experimen- 
tally observed "thermal masses" differ from the "band" 
masses by a renormalization factor of (1 + A), where A 
is the coupling constant for the interaction of electrons 
with phonons or other low-energy excitations. For Table 
I we used the values of A computed in the following way 
(see, e.g., Ref. 19): we assumed that the matrix elements 
of the electron-phonon interaction are constant within 
each of the 4 bands (a good approximation, see RefJ|), 
but different among the bands and for different interband 
transitions. If the matrix of the electron-phonon interac- 
tion is Uij, where i,j are the band indices, then the mass 
renormalization in the band i is 
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where Nj is the partial density of states per spin for the 
i-th band. Recall that the conventional Eliashberg cou- 
pling constant is A = UijNiNj / J2i ^i- The matrix 
U and the vector N calculated in Ref{j] were used to 
compute the fifth column in Table I. 

The agreement between the calculated and measured 
thermal masses can be characterized as excellent. Very 
importantly, this agreement is so good only because the 
calculated electron-phonon coupling differs by a factor of 
3 between the a and ir bands. This is the first direct 
demonstration of this important effect. The agreement 
between the calculated areas F and the experiment is 
also very good. Although Fi, F2 and F3 are overesti- 
mated by 35%, 15%, and 8%, respectively, the absolute 
values of these errors are only 0.5% (or less) of the total 
area of the corresponding Brillouin zone cross-sections. 
Even better appreciation of the significance of these er- 
rors can be gained from the observation that shifting the 
a band by 6.3 mRy down, and the tt* band by 5.5 mRy 
up brings the calculated areas to full agreement with the 
experiment. It is not at all clear whether or not such a 
small discrepancy with the experiment is meaningful. It 



is interesting, nevertheless, that after such an adjustment 
of the band positions the calculated masses agree with the 
experiment even better: for the three orbits in question 
the electron-phonon coupling constants deduced from the 
experiment by taking the ratio of the measured masses 
to the calculated masses are, respectively, 1.15,1.12, and 
0.43. After the Fermi level adjustment, they are 1.22, 
1.18, and 0.45. It is also worth noting that, for instance, 
a change in c/a ratio of 1.5% shifts the a and ir bands 
with respect to each other by 12mRy, or that a shift of 
the Fermi level by 6 mRy corresponds to a 0.05 e change 
in the number of electrons. This shows how sensitive the 
de Haas-van Alphen results are to the crystallography 
and stoichiometry. 

Another important observation reported in Ref.[l2| is 
the so-called "spin-zero". This is a suppression of the 
de Haas-van Alphen amplitude when the difference in 
the areas (in Tesla units) of the spin-split (by the ex- 
ternal field H) cross-sections is exactly H/2. This effect 
has been observed for orbit 8 in the field H = 17 T, 
when the field was tilted with respect to the crystallo- 
graphic axis by <f> — 15 — 18°. This means that (Fj — 
F 8 x )/cos(0) = 8.5 T, or AF 8 = F 8 T - Fg 1 « 8.1 T (note 
that the angle itself does not depend on the field in which 
the measurements are performed, but only on the Fermi 
surface geometry and Stoner renormalization) . It is easy 
to estimate this splitting in the first approximation, using 
the data from the Table I and the Stoner renormalization 
of 33%, calculated in Ref^l] A As = 2irmAE xc , where 
AE XC = 2/ibH(1 + S) is the induced spin-splitting of the 
bands near the Fermi level, enhanced by a Stoner factor 
(1 + S). This formula gives AF 8 « 7.1 T. A caveat here 
is that the induced spin-splitting need not be the same 
for all bands, in other words, while the average S is 0.33, 
individual «S"s may vary from orbit to orbit. To avoid 
this problem, we performed self-consistent LAPW calcu- 
lations in an external field of 1.8 kT (still well within the 
linear response regime) and measured dAAg,/dH explic- 
itly. Using these results, we found that for the actual 
field of 17 T AF 8 = 6.7 T, close to, but slightly smaller 
than the above estimate of 7.1 T. In other words, the cal- 
culated Stoner factor for this orbit is Sg = 0.26, smaller 
than the average over all bands, which is 0.33. Note 
that the experimental number of 8.1 T can be reconciled 
with the calculated mass, if S$ were sa 0.5, fairly close 
to the electron-phonon coupling constant for the same 
band, 0.47. We, however, believe that the coincidence is 
accidental, although we do not have any plausible expla- 
nation for the noticeable underestimation of the Stoner 
factor for this orbit. No "spin-zero" effect has been ob- 
served for the orbit 4, which has essentially the same mass 
as orbit 8. Our calculations for this orbit give AF4 = 6.9 
T; that is, the calculated Stoner factor for this orbit is 
S4 = 0.31. At the same time, the actual Stoner factor 
must be either larger than 0.60 or smaller than 0.18, for 
this orbit not to exhibit the "spin-zero" effect (this is ne- 
glecting deviations from a cylindrical shape, which are 
noticeably stronger expressed for this orbit than for the 
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TABLE I: Calculated de Haas-van Alphen parameters from present work (-Fcalc) compared to the experimental data (F cxp ) of 
RefJl| The masses are given in free electron mass units. 





Orbit 


-fcalc [T] 


m calc 


dm calc /dE [Ry" 1 ] 


A a 


|(l + A)m| 6 


F exp [T] 


|(l + A)m| cxp 


1 


a T-plane 


730 


-0.251 


1.1 


1.25 


0.56 


540 


0.54 


2 


a* T- plane 


1589 


-0.543 


2.7 


1.16 


1.17 






3 


7r T-plane 


34630 


1.96 


23 


0.43 


2.80 






4 


a A-plane 


1756 


-0.312 


1.2 


1.25 


0.70 


1530 


0.66 


5 


a* A-plane 


3393 


-0.618 


2.3 


1.16 


1.33 






6 


7r* T- plane 


31130 


-1.00 


4.1 


0.47 


1.47 






7 


7r TAM-plane 


458 


-0.246 


1.5 


0.43 


0.35 






8 


7r* TAM-plane 


2889 


0.315 


0.8 


0.47 


0.46 


2685 


0.45 



"Computed from Tables 1 and 2 of Ref.J 
'Computed from the preceding columns 



orbit 8). Further experimental studies on better samples 
should give more insight into this problem. 

Finally, we would like to discuss the problem of the 
"missing orbits". The amplitude of the de Haas- van 
Alphen signal is proportional toE3: 

_i X -chV^A irAF 



X = n 2 mc(l + \)k B T/heH, 

where I is the mean free path for the orbit in question. 
Thus, it is not surprising that the large orbits 3 and 6 are 
not observed; the Dingle exponent chV irA/eH£ is at least 
10 times larger than for the other orbits. However, the 
question remains for the orbits 2, 5, and 7. Let us start 
with the first two. We observe that, compared to the 
orbits 1 and 3, both Dingle factor and the thermal factor 
are reduced. The latter is smaller because the effective 
mass, m(l + A) is twice larger, which reduces the maximal 
temperature at which these orbits can be observed by 
a factor of two. The former is reduced because both 
the orbit size, y/A, is larger, and the mean free path, 
£ oc vf, (assuming the relaxation time is the same for 
both a and a* bands) is smaller (from Table 1 of Ref.|2l|, 
uf(ct) /vf{<J*) ~ 1-4). The total reduction of the Dingle 
exponent compared to orbit 4 is by a factor of 2 for orbit 
5, and of 1.4 for orbit 2. 

The absence of a signal from the orbit 7 seems puzzling. 
Its area and its thermal mass are the smallest of all or- 
bits, and the average velocity for this band is the highest 
(50% higher than for the a band). A very plausible ex- 
planation is that, as conjectured in Ref.^ and elaborated 
in Ref.^J, the impurity scattering rates differ drastically 
between the bands. If the dominant defects reside in 
the Mg plane (e.g., Mg vacancies), then such defects are 
very weak scatterers for the a bands for the simple reason 
that those bands have very little weight at the Mg atoms. 
However, this simple picture does not explain why orbit 
8, originating from the 7r* band, apparently has a small 
relaxation time and therefore is seen in experiment. Its 
velocity is close to (in fact, 15% smaller than) that of the 



7T band, its linear size is more than twice larger than that 
of orbit 7, so the scattering rate has to be at least 5 times 
larger. We do not have a plausible answer at the moment 
why the impurity scattering appears to be so suppressed 
for this orbit. Possibly, this is related to its parity (while 
the 7r band is even with respect to the z —> —z reflection, 
the 7T* band is odd). 

To conclude, we presented highly accurate calculations 
of the de Haas- van Alphen parameters for MgB2. Com- 
parison with the experiment reveals: (a) Absence of any 
mass (velocity) renormalization apart from that due to 
phonons. (b) A good agreement of the calculated cross- 
section areas with the experiment, (c) Excellent agree- 
ment of the calculated electron-phonon coupling with the 
dHvA mass renormalization, including very large dispar- 
ity between the coupling of the a— and 7r— bands, which 
clearly confirms the basic assumption of the two-gap 
model for superconductivity in MgB2. (d) Some under- 
estimation, despite a good qualitative agreement, of the 
calculated and measured Stoner factors for the it bands, 
(e) An indirect evidence of substantially different impu- 
rity scattering rates in the a and tt bands, (f) A problem 
which remains to be understood is the total suppression 
of the neck orbit, associated with the bonding tt band, 
given a clear observation of the much larger orbit from 
the electronically similar n* band. 

After this work was finished, we learned about similar 
works by Rosner et aLc9 and HarimaBa Their results, 
particularly those of Ref.|2^, are quite close to ours. Both 
paper employ similar methods and take full care of the 
k-mesh convergence. The remaining difference is a good 
gauge of how reliable are such calculations, in technical 
sense. 
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